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Abstract. Recently, the "(p+h)— wave" form of pairing symmetry has been proposed 
for the superconductivity in PrOs^Sb^ compound [Parker D, Maki K and Haas S, 
cond- mat/0407254]. In the present paper, a stationary Josephson junction as a 
weak-link between PrOs^Sb^ triplet superconductors is theoretically investigated. 
The quasiclassical Eilenberger equations are analytically solved. The spin and charge 
current-phase diagrams are plotted and the effect of misorientation between crystals 
on the spin current, and spontaneous and Josephson currents is studied. It is found 
that such experimental investigations of the current-phase diagrams can be used to 
test the pairing symmetry in the above-mentioned superconductors. It is shown that 
this apparatus can be applied as a polarizer for the spin current. 
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Transport of charge and spin in the weak link between misoriented PrOs^Sbu superconductor s2 
1. Introduction 

The pairing symmetry of the recently discovered superconductor compound PrOs^Sbyi 
is an interesting topic of research in the field of superconductivity [U 121 El- Supercon- 
ductivity in this compound was discovered in papers |H Ej and two different phases 
(A and B) have been considered for this kind of superconductor in Refs.jSJE!- Although 
authors of [2] at first considered the spin-singlet "(s + g)— wave" pairing symmetry for 
this superconductor, later it was specified that the spin-triplet is the real pairing sym- 
metry of the PrOs±Sb\2 complex P3E]. Using the Knight shift in NMR measurement 
authors of paper !8J, estimated the spin-triplet pairing symmetry for the superconductiv- 
ity in PrOs^Sbu- Consequently, the "(p+ h)— wave" model of the order parameter was 
proposed for the pairing symmetry of the superconductivity in PrOs^Sb^ compound, 
recently 0. In the paper |T], the self-consistent equation for the superconducting gap 
A (T) (BCS gap equation) has been solved for the finite temperature T numerically and 
for the temperatures T close to zero and the critical temperature T c analytically. For 
this compound, using the u (p+h)— wave" symmetry for the order parameter vector (gap 
function), the value of the A (0) has been obtained for both A and 5-phases, in terms 
of the critical temperature. In addition, the dependence of A (T) in the temperature 
limit of T — > and T — > T c have been obtained. Authors of paper pQ have investi- 
gated the temperature dependence of critical field, specific heat and heat conductivity. 
Also, the Josephson effect in the point contact between triplet superconductors with 
/—wave triplet pairing has been studied in Ref.jH]. In this paper the effect of misorien- 
tation on the charge transport has been studied and a spontaneous current tangential to 
the interface between the /-wave superconductors has been observed. Additionally, the 
spin-current in the weak-link between the /-wave superconductors has been investigated 
in our paper [TU]. In the paper [TU], this kind of weak-link device has been proposed 
as the filter for polarization of the spin-current. These weak-link structures have been 
used to demonstrate the order parameter symmetry in Ref. [TT] . 

In the present paper, the ballistic Josephson weak-link via an interface between two bulks 
of "(p+h)— wave" superconductor with different orientations of the crystallographic axes 
is investigated. It is shown that the spin and charge current-phase diagrams are totally 
different from the current-phase diagrams of the point-contacts between conventional (s- 
wave) superconductors > high T c (d-wave) superconductors and from the charge 
and spin-current phase diagrams in the weak-link between the /-wave superconductors 
EI]- We have found that in the weak-link structure between the "(p + h)— wave" 
superconductors, the spontaneous current parallel to the interface, as the characteristic 
of unconventional superconductivity, can be present. The effect of misorientation on 
the spontaneous, Josephson and spin currents for the different models of the paring 
symmetry (A— and B— phases in FigJTJ) are investigated. It is possible to find the value 
of the phase difference, at which the Josephson current is zero, but the spontaneous 
current tangential to the interface is present. In some configurations and at the zero 
phase difference, the Josephson current is not zero but has a finite value. This finite 
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value corresponds to a spontaneous phase difference, which is related to the misorien- 
tation between the gap vectors. Finally, it is observed that at certain values of the 
phase difference 0, at which the charge current is zero, the spin current is present and 
vise versa. In addition, in the configuration in which both gap vectors are directed 
along the c _L n axis (n is the normal to the interface unit vector), only the normal 
to the interface spin current s„ can be present and the other terms of the spin current 
are absent (s is the spin vector of electrons). Consequently, this structure can be used 
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Figure 1. A— phase (left), B— phase (right) order parameters and direction of a, b, c 
and z unit vectors (middle) 1 . A— and B— phases are high field (high temperature) 
and low field (low temperature) phases, respectively |14j . 



as a filter for polarization of the spin transport. Furthermore, our analytical and nu- 
merical calculations have shown that the misorientation is the origin of the spin current. 

The organization of the rest of this paper is as follows. In SecO we describe our 
configuration, which has been investigated. For a non-self-consistent model of the order 
parameter, the quasiclassical Eilenberger equations [T5] are solved and suitable Green 
functions have been obtained analytically. In Seed] the obtained formulas for the Green 
functions have been used for calculation the charge and spin current densities at the 
interface. An analysis of numerical results will be done in Sec0] The paper will be 
finished with some conclusions in Sec|SJ 



2. Formalism and Basic Equations 



We consider a model of a flat interface y = between two misoriented "(p + 
h)— wave" superconducting half-spaces (FigJ2J) as a ballistic Josephson junction. In the 
quasiclassical ballistic approach, in order to calculate the current, we use "transport- 
like" equations [T3] for the energy integrated Green matrix g (y F ,r, e m ) 

(1) 



e m a 3 + iA,g 



0. 



and the normalization condition gg = 1, where e m = nT(2m+l) are discrete Matsubara 
energies m — 0, 1, 2..., T is the temperature, vj? is the Fermi velocity and 03 = 03 ® I in 
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Figure 2. Scheme of a flat interface between two superconducting bulks, which are 
misoriented by angle a. In geometry (i), the afo-plane on the right side is rotated as 
much as a around the c-axis and in geometry (ii) the c-axis on the right side is rotated 
around the 6-axis. ab~ planes and c— axis have been shown in FigQ 

which &j (j = 1,2,3) are Pauli matrices. Also the matrix structure of the off-diagonal 
self energy A in the Nambu space is 

A=( ° d ^ 2 V (2) 

where, the gap vector d, is the three dimensional counterpart of the BCS energy-gap 
function A for the case of triplet superconductivity. The Green matrix g can be written 
in the form: 

" _ [ 9l + (92 + g2^) io-2 

\ i&2 (93 + gs<r) 1*2(91 ~ gio-)icr 2 

Here, ^i+gi5" and (^2+g20")^2 [i*2(gz+^a)] are normal and anomalous Green function 
matrix, respectively. The terms of gi and are coefficients of Green function matrix 
expansion in terms of the unit matrix and Pauli 2x2 matrices in the Nmabu space. 
Also, the terms g\ and gi determine the charge and spin current densities through the 
equations (jHJ) and (jSj), respectively and g 2 and g 3 are used to determine the gap vector, 
d, using the self consistent relations (Eq|S]and it's conjugate). It is remarkable that, in 
this paper, the unitary states, for which d x d* = 0, is investigated. Also, the unitary 
states vectors d 12 can be written as 

d n = A„exp^ n , (4) 

where Ai^ are the real vectors in the left and right sides of the junction. The gap (order 
parameter) vector d has to be determined from the self-consistency equation: 

d (v F , r) = 2vrTiV (0) ^ (V (v F , v^) g 2 (v^, r, e m )) (5) 

rn 

where V (v>, v^), is a potential of pairing interaction, (...) stands for averaging over the 
directions of an electron momentum on the Fermi surface v' F and N (0) is the electron 
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Figure 3. Component of charge current normal to the interface (Josephson current) 
versus the phase difference <f> for A and B— phases, T/T c — 0.08 and different 
misorientations. Part (a) is plotted for geometry (i) and part (b) for geometry (ii). 
Charge currents are given in units of jo = § c/V(0)ufAo(0). 



density of states at the Fermi level of energy. Solutions to Eqs. (HJ) and © must satisfy 
the conditions for Green functions and vector d in the bulks of the superconductors far 
from the interface as follow: 

- I _L_ \ £m °~ 3 + ^2,1 (R\ 

g{±oo) = (6) 

d (±oo) = A 2 ,i (v F ) exp (t 1 ^ + #2,1 J , (7) 

where <fi is the external phase difference between the order parameters of the bulks. 
Eqs. (0) and (0) have to be supplemented by the continuity conditions at the interface 
between superconductors. For all quasiparticle trajectories, the Green functions satisfy 
the boundary conditions both in the right and left bulks as well as at the interface. 
The system of equations Qand (0) can be solved only numerically. For unconventional 
superconductors such solution requires the information of the function V(vp,v^). 
This information, as that of the nature of unconventional superconductivity in novel 
compounds, in most cases is unknown. Usually, the spatial variation of the gap 
vector and its dependence on the momentum direction can be separated in the form 
of A(vp, y) = A(vp)\I/(y). It has been shown that the absolute value of a self-consistent 
order parameter and ^(y) are suppressed near the interface and at the distances of the 
order of the coherence length, while its dependence on the direction in the momentum 
space (A(vp)) remains unaltered [TH] . Consequently, this suppression doesn't influence 
the Josephson effect drastically. This suppression of the order parameter keeps the 
current-phase dependence unchanged but, it changes the amplitude value of the current. 
For example, it has been verified in Refs. [T5] for the junction between unconventional 
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d-w&ve, in Ref . [16 for the case of "/-wave" superconductors and in Ref . ^7] for pinholes 
in 3 He that, there is a good qualitative agreement between self-consistent and non-self- 
consistent results. Also, it has been observed that the results of the non-self-consistent 
investigation of D—N—D structure in JH] are coincident with the experimental results of 
[TH] and the results of the non-self-consistent model in j20j are similar to the experiment 
|21j . Consequently, despite the fact that self-consistent numerical results cannot be 
applied directly for a quantitative analysis of the real experiment, only a qualitative 
comparison of calculated and experimental current-phase relations is possible. In our 
calculations, a simple model of the constant order parameter up to the interface is 
considered and the pair breaking and the scattering on the interface are ignored. We 
believe that under these strong assumptions our results describe the real situation 
qualitatively. In the framework of such model, the analytical expressions for the current 
can be obtained for an arbitrary form of the order parameter. 



3. Analytical results. 

The solution of Eqs.(HJ) and (j3J) allows us to calculate the charge and spin current 
densities. The expression for the charge current is: 

j e (r) = 2meTN (0) ^ (v^i (v F , r, e m )> , (8) 



and for the spin current we have: 



j Si (r) = 2ittAtN (0) V (v F (e igl (v F , r, e m ))) (9) 



where, e$= (x, y, z) . We assume that the order parameter does not depend on the 
coordinates and in each half-space it equals its value (J7J) far from the interface in the 
left or right bulks. For such a model, the current-phase dependence of a Josephson 
junction can be calculated analytically. It enables us to analyze the main features of 
current-phase dependence for different models of the order parameter of u (p + h)— wave" 
superconductivity. The Eilenberger equations (JJ) for Green functions g, which are 
supplemented by the condition of continuity of solutions across the interface, y — 0, and 
the boundary conditions at the bulks, are solved for a non-self-consistent model of the 
order parameter analytically. Two diagonal terms of the Green matrix which determine 
the current densities at the interface, y = 0, are shown below. For the term relating to 
the charge current we obtain: 

= eJSh + ^2) cos/3 + ^(^1^2 + £m) sin/3 
91 U i7/e m (fii + fla) sin/5 + (fii0 2 + el) cos/3 + AiA 2 ' 1 ' 

and for the case of spin current we have: 

gl (0) = Ax x A 2 

(B - l) 2 (rA + e m )(r}fl 2 + e m ) exp(z/3) - (B + 1) 2 (t]Q 2 - e m )(»A - £ m) exp(-i/3) 



2r ? (A + J B)|A 1 | 2 |A 2 | 2 



(11) 
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Figure 4. Components of charge current tangential to the interface versus the 
phase difference cf> for A and B— phases, geometry (ii), T/T c = 0.08 and the different 
misorientations. Part (a) is plotted for a;— component and part (b) for z— component. 



where r] = sgn (v y ), fi n = \J e 2 m + | A n | , j3 = ip x - ip 2 + (f>, 

B _ Tie m (Vti + Q 2 ) cos (3 + ijp^h + g§J sin/3 
~ ir)£ m {&i + fi 2 ) sin/3 + + e 2 m ) cos (5 + A X A 2 



and 



(B - 1) exp{i/3) (B + 1) exp(-i/3) 



(13) 

(^1 - £ m )(^ 2 - £m) ' fafil + ^m)(^2 + £m)J ' 

Also, n = 1,2, label the left and right half-spaces respectively. We consider a rotation 
R only in the right superconductor (see, Fig|2J), i.e., d 2 (k) = _Rdi(^ _1 k); k is the 
unit vector in the momentum space. The crystallographic c-axis in the left half-space 
is selected parallel to the partition between the superconductors (along the z-axis in 
Figs ll|2|) . To illustrate the results obtained by computing the formula (fT0|) . we plot the 
current-phase diagrams for different models of the "(p + h)— wave" pairing symmetry 
(114115)1 and for two different geometries. These geometries correspond to the different 
orientations of the crystals on the right and left sides of the interface (FigEJ: 

(i) The basal afe-plane on the right side has been rotated around the c-axis by a; Ci||c2. 

(ii) The c-axis on the right side is rotated around the 6-axis by a (y-axis in FigEJ; 
bi||b 2 . 

Further calculations require a certain model of gap vector (order parameter vector) d. 



4. Analysis of numerical results 

In the present paper, two forms of "(p + h)— wave" unitary gap vector d in PrOs^Sb^ 
are considered. The first model to explain the properties of the A-phase of PrOs^Sb^ 
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Figure 5. Components of spin current (s y ) tangential to the interface versus the phase 
difference (f> for geometry (ii), Sr = 0.08 and different misorientations between the A 
and B— phase of "(p + h)— wave." Part (a) is plotted for component and part (b) 
for 2— component. Spin currents are given in units of j S Q — j1iN(0)vfAq(0). 



is (left side of FigflJ: 

d = A (T)(k x + ik y )^{l -k\-k\- kj)z. (14) 

The coordinate axes x, y, z are chosen along the crystallographic axes a, b, c on the left 
side of Fig |21 k is the unit vector along vp. The scalar function A = A (T) describes 
the dependence of the gap vector d on the temperature T. The second model to describe 
the gap vector of the -B-phase of PrOs^Sb^ is (right side of FigJTJ): 

d = A (T)(k x + tk y )(l-kt)z. (15) 

Our numerical calculations are done at the low temperatures, T/T c = 0.08, and 
we have used the formulas ln(A(T)/A(0)) = — 7n ^ (^oy) 3 f° r the A— phase and 
ln(A(T)/A(0)) = - 13 S 3l (^o)) 3 for the 5-phase, from the paper [H], for temperature 
dependence of the gap functions A (T) at the low temperatures (T << T c ). Also, in 
the paper j3] the value of A (T) has been calculated in terms of the critical temperature 
for both A and B- phases. They are A(0)/T c = 2.34 and A(0)/T c = 1.93 for A and 
B— phases, respectively. Using these two models of order parameters f!14| ITH|l and 
solution to the Eilenberger equations (fTUl l3~j). we have calculated the current density at 
the interface numerically. These numerical results are listed below: 
1) In part (a) of FigJBl the component of current normal to the interface of current, which 
is known as the Josephson current, is plotted for both A and B— phases, geometry (i), 
misorientations a = it/4 and a = tt/6. It is observed that the critical values of current 
for the B— phase is larger than the A— phase. Also, unlike Josephson junction between 
the conventional superconductors, here, at = 0, the current is not zero. The current is 
zero at the phase difference value <fi = <po, which depends on the misorientation between 



Transport of charge and spin in the weak link between misoriented PrOs^Sb\2 superconductors^ 




Figure 6. Component of charge and spin currents (J y and j sy ) normal to the interface 
versus the phase difference <j) for geometry (ii), rf- = 0.08, a = -| and a = \. Part (a) 
is plotted for A— phase and part (b) for B— phase. 



the gap vectors. In FigJ3J the value of the spontaneous phase difference 0o is close to 
misorientation a. 

2) In part (b) of FigOU the Josephson current is plotted for both A and B— phases, 
geometry (ii) and different misorientations. Again, the maximum value of the current 
for the B— phase is larger than for A— phase. Increasing the misorientation between the 
gap vectors, the maximum value of the current decreases. It is demonstrated that at the 
phase difference values <fi = 0, <fi = it and <ft = 2ir, the Josephson current is zero while 
both spontaneous and spin currents are not zero and have a finite value. Increasing the 
misorientation between the gap vectors decreases the derivative of the current with 
respect to the phase difference close to = 7r, which is known as the SQUID sensitivity 
and it is important from the application point of view. 

3) In Fig0J the tangential components of the charge current (x and z— components) in 
terms of the phase difference <fi are plotted. It is seen that at cf) = 0, <p — ft an d <p = 2vr, 
at which the Josephson current is zero, the parallel spontaneous currents have finite 
values. Although the normal component of charge current (see part (b) of FigOJ) is an 
odd function of the phase difference with respect to the line of = it while the parallel 
charge currents for this geometry (ii) are even functions of the phase difference with 
respect to <fi = tc (compare part (b) of Figl^Jwith Fig0J). 

4) In FigEl the tangential components of the spin (s y ) current are plotted in terms of 
the phase difference, for geometry (ii) and different misorientations. By increasing the 
misorientation the maximum value of the spin current increases. In spite of the charge 
current for this state, the spin current at the phase differences <fi = 0, <p — tt and <fi = 2ir 
is zero exactly(compare Fig0]with FigEI) • 

5) In FigJBl the normal component of charge and spin current (j y and j s ) are plotted 
for different misorientations and A and B— phases respectively. An interesting case in 
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our observations, is the finite value of normal spin current at = 0, = 7r and = 2n 
at which the normal charge current (j y ) is zero (see FigJSJ). 

6) In part (a) of FigEl the Josephson current is plotted in terms of the phase difference 
for the case of p— wave, and A and B— phases of "(p + h)— wave". The p— wave pairing 
symmetry as the first candidate for the superconducting state in Sr 2 Ru04 is as follows 
J231: 

d = A (T)(k x + ik v )z (16) 

It is observed that the maximum value of the Josephson current (j y ) of the junction be- 
tween the p— wave superconductors, is larger than for the B— phase of "(p + h)— wave" 
and the Josephson current of B— phase is larger than its value for the A— phase coun- 
terpart. Also, the place of the zero of the charge current for these three types of 
superconductor (geometry (i)) is the same. It is at the spontaneous phase difference 
which is close to the misorientation (j) = a (see, part (a) of FigOand part (a) of FigEJ). 

7) In part (b) of FigJTJ the normal component of the spin current is plotted for p— wave, 
A and B— phases of "(p + h)— wave" pairing symmetries and for a specified value of 
misorientation a — f. In both (a) and (b) parts of FigJTJ the maximum value of the 
current of junction between the p— wave superconductors is larger than for the B— phase. 
For B— phase the current has the maximum value, which is larger than the value for 
junction between the "(p + h)— wave" superconductors in the A— phase. This different 
characteristic of the current-phase diagrams enables us to distinguish between the three 
states. Also, it is observed that at phase differences = 0, = it and = 27r, the spin 
current has a finite value which may have its maximum value. This is a counterpart of 
part (b) of FigEJ in which the charge currents are zero at the mentioned values of the 
phase difference, but the spin current has a finite value. 

Furthermore, our analytical and numerical calculations have shown that the origin 
of the spin current is misorientation between the gap vectors [cross product in Eq.flHJ)]. 
Because the geometry (i) is a rotation by a around the z— axis and both of the left 
and right gap vectors are in the same direction, the cross product between gap vectors 
and consequently, the spin current for geometry (i) is zero. It is shown that in this 
structure (y— direction is normal to the interface, c— axis is selected in the z— direction 
and rotation is done around the y— direction) only the current of the s y flows and other 
terms of the spin current are absent. So, this kind of weak-link experiment can be used 
as a filter for the polarization of spin transport. Since the spin is a vector, the spin 
current is a tensor and we have the current of spin s y in the three x, y and z directions. 

5. Conclusions 

We have theoretically studied the the spin and charge transport in the ballistic Josephson 
junction in the model of an ideal transparent interface between two misoriented 
PrOs±Sbi2 crystals with u (p+h)— wave" pairing symmetry, which are subject to a phase 
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Figure 7. Component of charge and spin currents normal to the interface versus 
the phase difference (f>, S- = 0.08, a = ?, A— phase, B— phase and p— wave pairing 
symmetry. Part (a) is plotted for charge current and geometry (i) and part (b) is for 
the case of spin current and geometry (ii) 



difference 0. Our analysis has shown that the different misorientations and different 
models of the gap vectors influence the spin and charge currents. This has been shown for 
the charge current in the point contact between two bulks of "/—wave" superconductors 
in and for the spin current in the weak link between "/-wave" superconductors in 
[TP] . In this paper, it is shown that the misorientation of the superconductors leads to a 
spontaneous phase difference that corresponds to the zero Josephson current and to the 
minimum of the weak link energy. This phase difference depends on the misorientation 
angle. We have found a spontaneous charge current tangential to the interface which 
is not equal to zero in the absence of the Josephson current generally. It has been 
found that the spin current is the result of the misorientation between the gap vectors. 
Furthermore, it is observed that a certain model of the gap vectors and geometries can 
be applied to polarize the spin transport. Finally, as an interesting and new result, it 
is observed that at certain values of the phase difference 0, the charge-current vanishes 
while the spin-current flows, although the carriers of both spin and charge are the same 
(electrons). The spatial variation of the phase of the order parameter plays a role as 
the origin of the charge current and, similarly, due to the broken Q s P m -° rblt symmetry, a 
spatial difference of the gap vectors in two half-spaces is the cause of spin currents. This 
is because there is a position-dependent phase difference between "spin up" and "spin 
down" Cooper pairs and, although the total charge current vanishes, there can be a net 
transfer of the spin. Therefore, in our system, there is a discontinuous jump between 
the gap vectors and, consequently the spin currents should generally be present. For 
instance, if up-spin states and down-spin states have a velocity in the opposite direction, 
the charge currents cancel each other whereas the spin current is being transported. 
Mathematically speaking, j c har g e = j| +j|?js P in = j| — j 4. , so it is possible to find the 
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state in which one of these current terms is zero and the other term has a finite value 
[22] • in conclusion, the spin current in the absence of the charge current can be observed 
and vice versa. 
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